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ABSTRACT: Previous large N calculations are combined with numerical work at 
N = 4 to show that the Minimal Standard Model will describe physics to an accuracy 
of a few percent up to energies of the order 2 to 4 times the Higgs mass, Mjj, only if 
M^i — 710 ± 60 GeV. This bound is the result of a systematic search in the space of 
dimension six operators and is expected to hold in the continuum. Given that studying 
the scalar sector in isolation is already an approximation, we believe that our result is 
sufficiently accurate and that further refinements would be of progressively diminishing 
interest to particle physics. 



1. INTRODUCTION AND CONCLUSION. 



Our goal is to obtain an estimate for the triviality bound on the Higgs mass in the mini- 
mal standard model. Much work has preceded this paper (for examples consult [1,2,3,4,5,6] 
and the review [7]). We build on these results and generalize them. Specifically, we deal 
more systematically with the quantitative uncertainty resulting from the arbitrariness of 
the coefficients of higher dimensional operators in the scalar sector [8] . By doing this we 
aim to obtain a number that is meaningful beyond lattice field theory and directly relevant 
to particle physics. 

The overall framework of the approach has been reviewed before (e.g. [8]) and will not 
be repeated here. Its main simplifying feature is to treat the scalar sector of the minimal 
standard model in isolation. When considering further efforts in this framework the po- 
tential impact on particle physics should be evaluated against the accuracy of neglecting 
other interactions, e.g. with the top quark. We believe that our result is reliable to a 
reasonable degree, given that the whole framework is an approximation, and we feel that 
further refinements would be of progressively diminishing interest to particle physics*. 

Our result is that the minimal standard model will describe physics to an accuracy 
of a few percent up to energies of the order 2 to 4 times the Higgs mass, Mjj, only if 
Mjj < 710 ± 60 GeV. The two major assumptions made are that ignoring all couplings 
but the scalar self-coupling is a good approximation and that any higher energy theory 
into which the standard model is embedded will not conspire to eliminate all dimension 
six scalar field operators in the low energy effective action. The number we obtain is not 
surprising because of its closeness to tree level bounds [1,3]; what has been achieved is to 
finally show that in any reasonable situation higher orders in perturbation theory cannot 
change it substantially although quite strong scalar self-interactions are possible. 

In the sequel we shall rely quite heavily on [9] but we also try to make the paper 
accessible to readers who are not familiar with [9]. Our basic strategy was to first use 
expansion techniques in the generalization of the 0(4) symmetric scalar sector to O(N) 
to obtain an analytical non-perturbative estimate for the bound and then follow up with 
Monte Carlo simulations at the physical value N = A. In the next section we summarize 
needed information at N = oo and add some new results in this limit. The following 
section presents our numerical work with some emphasis on the checks that were made to 
ascertain control over systematic errors. The last section explains our main conclusion. 
Appendix A gives a few more details on the new large N results, and finally in appendix 
B we collect the numbers obtained from our simulations in several tables matching the 
graphs shown in the main text. 

* An exception would be further investigations of the Higgs width on the lattice. 
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Our notational conventions are: The Higgs mass in physical units (GeV), defined as 
the real part of the resonance pole, is denoted by Mjj and the width by T^. The same 
quantities in lattice units are denoted by lower case letters, mjj and 7^. The matrix 
element of conventionally normalized currents of broken symmetries between the vacuum 
and single Goldstone boson states (referred to as pions, n) is denoted by F in physical 
units and by / in lattice units. The scalar self-coupling is defined by: 

9 = ^f ■ (1-1) 

In the N = 00 section everything is written in terms of the above N = 4 notation. We use 
A to denote a generic cutoff. 

Most of our work is on the F4 lattice which can be thought of as embedded in a 
hypercubic lattice from which odd sites (i.e. sites whose integer coordinates add up to an 
odd sum) have been removed. The lattice spacing of the hypercubic lattice is a. It is set 
to unity when mjj, 7^ and / are used. The F4 lattice is always fully symmetric having, 
when finite, L sites in each principal axis direction so that the total number of sites is L 4 . 
Usually, x,x',x" denote sites, < x,x' >,/,/' links, « » next-nearest-neighboring 
pairs, </,/'> pairs of links, and the field is constrained by $ 2 (x) = 1. 



2. RESULTS AT LARGE N. 



Here we summarize the large N results of [9] relevant to our numerical work. They 
contain predictions that can be directly compared to AT = 4 Monte Carlo data and eval- 
uations of observable cutoff effects on n — n scattering. The dependence of the bound 
on the magnitude of the observable cutoff effects is relatively insensitive: a change by a 
factor of 3 induces a variation of 50 GeV in the bound in the worst case. Thus, we do not 
worry about jj- corrections to the cutoff effects. Cutoff effects could also be calculated in 
perturbation theory but we have argued (see Appendix B of [9]) that the computation 
is probably more reliable. Hence we use the results here. 

2.1. Relaxing the bound. 

The naive expectation that heavier Higgs masses are obtained when the bare scalar 
self-coupling is increased is upheld by nonperturbative calculations. The search for the 
bound can therefore be restricted to nonlinear actions. 

Among the nonlinear actions the bound is further increased by reducing as much as 
possible the attraction between low momentum pions in the I = J = channel. Given 



3 



a bare action the approximate combination of parameters achieving this is identified as 
follows. Expand in slowly varying fields and use a field redefinition to bring the action, 
including terms up to fourth order in the momenta, to the form 




2N 



2N 



{d^-d v $- l -5^ v dj-dj) 2 . (2.1.1) 



At N = oo 62 has no effect and the bound depends monotonically on b\, increasing with 
decreasing b\. Overall stability of the homogeneous broken phase restricts the range of b\. 
For example on the F4 lattice, at the optimal value of b\ the bound is increased by about 
100 GeV relative to the simplest non linear action. 

The rule in the above paragraph does not lead to an exactly universal bound. Different 
bare actions that give the same effective parameter b\ can give somewhat different bounds 
because the dependence of physical observables on the bare action is highly nonlinear. 
For example, at the optimal b\ value, Pauli-Villars regularizations lead to bounds higher 
by about 100 GeV than some lattice regularizations. This difference between the lattice 
and Pauli-Villars can be traced to the way the free massless inverse Euclidean propagator 
departs from the 0(p 2 ) behavior at low momenta. For Pauli-Villars it bends upwards to 
enforce the needed suppression of higher modes in the functional integral, while on the 
lattice it typically bends downwards to reflect the eventual compactification of momentum 
space. 

Because we desire to preserve Lorentz invariance to order 1/A 2 we use the F4 lattice 
and, on the basis of the above observations, there are three stages of investigation. The 
first stage is to investigate the naive nearest-neighbor model. This should be viewed as the 
generic lattice case where no special effort to increase the bound is made. Since this case has 
been investigated thoroughly in [5] we can proceed to more complicated actions with well 
tested methods of analysis. The next stage is to write down the simplest action that has a 
tunable parameter b\ . The last stage is to add a term to eliminate the "wrong sign" order 
term in the free nearest-neighbor propagator, amounting to Symanzik improvement of 
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the large N pion propagator. The three F4 actions, investigated at large N, are given by 
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These actions were chosen because they have relatively simple large N limits and permit 
a simultaneous study of both F4 and hypercubic lattices. In each case, at constant 
Pq is varied tracing out a line in parameter space approaching a critical point from the 
broken phase. This line can also be parameterized by mjj or g. For actions S'2 and S'^ 
is chosen so that on this line the bound on Mjj is expected to be largest. A simulation 
produces a graph showing ^j- as a function of rrifj along this line. The y-axis is turned 
into an axis for Mjj by Mjj = 1 ^j L x 246 GeV. The large predictions for these graphs 
are shown in Figure 2.1. Since action S3 was not treated in [9] we include a brief account 
in appendix A. 



2.2. Cutoff Effects. 



At N = 00 the cutoff effects to order 1/A 2 on the Higgs width and n — n scattering 
are parameterizable by: 

Seff = S R (g) +cexpi-^ 2 /9Mg) ■ (2.2.1) 

Sr contains only universal information and so does O. All the non-universal information 
is in the g independent parameter c. The form of S e jj reflects the factorization of cutoff 
effects into a universal g dependent function and a g independent non-universal amplitude. 
The bound is increased by first varying the non-universal part so that c decreases at 
constant g and then going along the selected line to higher g. 

5|^4|2 denotes the fractional deviation of the square of the n a — n a (a identifies one of the 
N — 1 directions transverse to the order parameter in internal space) scattering amplitude 
at 90 degrees in the CM frame at energy W from its large N universal value. A plot of 8^2 
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Figure 2.1 The Higgs mass Mjj = ^j 1 x 246 GeV in physical units vs. the Higgs mass 
ran in lattice units for the three actions, eq. (2.1.2). 

as a function of Mjj{mjj) for the three actions is presented in Figure 2.2. If one considers 
only the magnitude of cutoff effects as a function of mjj = one might conclude that 
the bound obtained with action would be larger than the bound obtained with S^. This 
conclusion proves to be wrong when the mass in physical units is considered. The values of 
Mjj in GeV, determined from Mjj = ^P- x 246 GeV, are put on three horizontal lines in 
Figure 2.2 at 8\M2 = 0.005, 0.01, 0.02. At large iV the bound increases when going from 
S[ to ^ and then to S3 by a little over 10% at each step. For example, for 5\a\2 = .01 we 
get bounds on Mjj of 680, 764, 863 GeV for S[, and £3 respectively. 




2.3. Width. 



The width T# is important for phenomenology and for lattice work (see sect. 3.4). 
With massless pions perturbation theory seriously underestimates T# when Mjj is large 
[9]. 
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Figure 2.2 Leading order cutoff effects in the invariant n — n scattering amplitude at 
90° at center of mass energy W = 2Mjj vs. the Higgs mass in lattice 
units for the three actions. The values of Mjj in GeV determined from 
Mh = —f 1 x 246 GeV are put on the three horizontal lines at 8\a\2 = 
0.005, 0.01, 0.02. 



It would therefore be desirable to determine the width non-perturbatively. Up to date 
the only methods known for measuring the width within a numerical simulation require a 
non-zero pion mass [10]. To study the effects of a non-zero pion mass we computed the 
width in this case at large N. The mass for the pions was induced by an external magnetic 
field that breaks the symmetry explicitly [11]. Because the cutoff effects on the width are 
very small we show in Figure 2.3 only the universal part and compare it to the leading 
order perturbative values. Shown are the two results for M n = 0, 100, 200, and 300 GeV. 

We see that the deviations from perturbation theory decrease when the pion mass 
increases. Therefore to detect nonperturbative effects on the width at N = 4 one would 
have to deal either with the massless case directly or, at least, with quite light pions, say 
M n < j*Mfj. This is one point we believe deserves further study. 
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Figure 2.3 The universal part of the width vs. M^j. The solid line displays the large N 
result scaled to N = 4 and the dotted line shows the leading order term in 
perturbation theory. From left to right the lines correspond to pion masses 
Mjr = 0, 100, 200, and 300 GeV. 

3. THE PHYSICAL CASE A = 4. 

The primary aim of the numerical work is to produce at A = 4 the analogues of the 
graphs in Figure 2.1. The actions simulated are slightly different than those investigated 
at large A. Firstly, the factors A accompanying the couplings in (2.1.2) are omitted. 
Secondly, we take advantage of the fact that on the F4 lattice, unlike on the hypercubic 
lattice, the simplest action with a tunable parameter b\ can be constructed in a way that 
maintains the nearest-neighbor character of the action. This is done by coupling fields 
sited at the vertices of elementary bond-triangles. Finally we add a term to eliminate 
the "wrong sign" order term in the free nearest-neighbor propagator, leaving, unlike 
in S'^, the term with the tunable b\ coupling unchanged. This new term couples next- 
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nearest-neighbors and amounts to tree level Symanzik improvement.* The three F4 actions 
simulated (with S{ having the same expansion in slowly varying fields as S[ investigated 
at large N) are 

S 1 = - 2/3 ®( x ) ■ *(*0 

<x,x'> 

($(x) ■ <S(a/)) ■ 

x <W> 

1,1'nxjt®, inx'jtQ, I'nx"^® 

x,x f ,x ff all n.n. ^ 

S 3 = -2(2/3 +/3 2 ) J2 + (Po + fo) E 

<x,x'> <<Cx,x'~^> 

X <ll>> 

i,i'nx^(ti, inx'=i<», l'nx"ji$ 

x.x r .x ff all n.n. 

The data for action S\ can be found in [5] . Preliminary data for action £2 were presented 
in [8,12] and references therein. Our main task is to finalize the results for action S2 and 
carry out some new measurements for action S3. 

Each action is investigated in two main steps. First, the phase diagram is established; 
next, a particular line is chosen in the broken phase, which for actions £2 an d S3, amounts 
to picking a value for /?2- On this line we make several measurements at different values of 
flo approaching the critical point Podfo)- We chose P2 based on the large iV criteria and 
on the measured structure of the phase diagram. We may therefore be missing the "best" 
action, but, from our experience at large N, we expect at most a 20 — 30 GeV additional 
increase in the bound. Because we object to excessive fine-tuning and since one may view 
what we are doing already as some amount of fine tuning, the 20 — 30 GeV might go in 
either direction and represents the systematic uncertainty that we assign to the question 
of fine tuning. 

3.1. Methods in General. 

We follow closely the approach of [5]. We use a Metropolis algorithm to map out the 
phase diagram and a single cluster spin reflection algorithm, tested against the Metropolis 
algorithm, for the actual measurements. Typically we use 10, 000 — 100, 000 lattice passes, 
depending on lattice size and couplings, and simulate systems of increasing sizes with even 

On the hypercubic lattice tree level improvement of the simplest action may also help 
reduce Lorentz violation effects at order 1/ A 2 and was investigated in [6] . 
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L (to avoid frustration effects), L = 6, • • • , 16. The spin update speeds on a single processor 
CRAY Y-MP are about 20^ for action Si, 50^ for S 2 and 55^f for S3. 

The total computer time invested is 400 hrs for action S\, 500 hrs for £2 and 400 hrs 
for S3. The approximate constancy of the total amount of time spent for each action, in 
spite of the significant increase in complexity indicated by the rise of time per spin update, 
reflects the diminishing need for checks as confidence in the numerical methods builds up. 
The total amount of time spent, approximately 1300 hrs, is quite modest and shows that 
careful preparation, continuous support by analytical work and an incremental approach 
pay off. 

Statistical errors are always treated by the Jackknife method and in least x 2 -fits cor- 
relations between measurements are taken into account. 

3.2. Phase Diagrams. 

For actions £2 and S3 we determine critical points fa c {fa) a t several fixed fa's using 

y ®(x) 

Binder cumulants of the magnetization M. = L 4 — • The critical points are obtained 
from the intersection points of the Binder cumulant graphs for different volumes, produced 
by reweighting and patching [13] the measurements from several couplings fa. The results 
are shown in Figures 3.1 and 3.2 and the actual numbers are listed in Tables 1 and 2. 
They compare well with our large N results which were helpful in deciding where to scan 
in the first place. At large N the critical lines are straight, a feature that seems to persist 
at N = 4. 

Based on our large N work we know that we want extremal values of fa [9]. By 
this we mean the most negative value of fa for which the leading term in the long wave- 
length expansion, eq. (2.1.1), still has the usual ferromagnetic sign on the critical line (i.e. 
fac + fa <i 0). We therefore select fa = —0.11 for actions S2 and S3. The critical points 
at fa = —0.11 are fa c = 0.1118(3) and fa c = 0.1113(2) respectively. The specific points 
we choose to study in detail are shown in Figures 3.1 and 3.2. 

3.3. Coupling Constant. 

To obtain the coupling constant g, eq. (1.1), we need / and mjj. f is obtained from 
measuring the magnetization, M., and the pion wave function renormalization constant, 
Z^, defined from the residue of the pole at zero momentum in the pion-pion propagator. 
The pion field, 7?, is defined as the component of $ transverse to Ai, and for low momenta 
we have < |7?(p)| >= -^f + regular terms. The magnetization is obtained by extrapolating 
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Figure 3.1 Phase diagram for action $2 • The solid line is a least square straight line fit 
to the data that are denoted by diamonds. The squares indicate the points 
where we made simulations to determine ran an d fir- They all lie on the 
vertical line = —0.11. 

to L = 00 the quantities < Ai^ >l with an 0(1/ L 2 ) correction, using the methods in [14]. 
/ is then obtained from / = Ai/Z n . This method is safe and the reliability of the numbers 
one obtains for / is high for our purposes. The error never exceeds 1.5% in the region of 
higher Higgs masses which we are interested in. The estimate of / by analytical methods 
[3] and [15] has an error of order 5% in the same region; thus / is better determined by 
Monte Carlo. This is due mainly to the good theoretical control one has over finite volume 
effects [5,6,14]. Unfortunately, for the determination of ran we are not so lucky. 

To obtain ran we measure the correlations of zero total three-momentum sigma states 
and low relative momentum two-pion states at different time separations. The sigma 

— * 

field, a, is defined as the component of $ parallel to M. ran is the energy of one of the 
lightest states created by superpositions of the a field and two-pion composite operators 
at zero total three-momentum and is obtained from the eigenvalues of the measured time 
correlation matrix. The evaluation of ran is n °t as clean as that of / and the next 
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Figure 3.2 Same as Figure 3.1 but for action S3. 

subsection discusses the determination of the numbers in greater detail. 

The main result of this paper is in Figure 3.3 which shows Mjj = 246a/ g/S GeV as 
a function of ni}j for all three actions. The actual numbers are given in Tables 3, 4 and 
5. One clearly sees the progressive increase of the bound. A glance at Figure 2.2 shows 
that in all cases the cutoff effects on the pion-pion scattering are below a few percent 
even at the maximal Mjj of each curve. Thus we can take the largest of these maxima as 
our bound. The ordering of the points and their relative positions are in agreement with 
Figure 2.1, while the differences in overall scale, reflecting the difference between N = 00 
and N = 4, come out as expected [9]. 



3.4. More about mass determination. 



The main source of systematic errors is in the evaluation of rriff. In an infinite volume 
the Higgs particle would decay into two pions and extracting m# from the fall-off of the 
a-a correlation function would yield nonsense. In a finite volume the decay is prohibited 
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Figure 3.3 The Higgs mass Mjj = ^j- x 246 GeV in physical units vs. the Higgs 
mass mjj in lattice units from the numerical simulations. The diamonds 
correspond to action S\ [5], the squares to action S2 and the crosses to 
action S3. 

or severely restricted by the rather large minimal amount of energy even the softest pions 
have due to momentum quantization. This makes the measurement possible but leaves 
one with the difficult task of estimating the accuracy of the so determined real part of the 
resonance pole. 

In [5] we dealt with action Si. As a check of the results obtained in the broken phase 
we also measured the coupling g in the symmetric phase (as defined from the zero four- 
momentum four-point correlation) and used perturbation theory to predict g in the broken 
phase. This method is free of finite width contamination and gave results consistent with 
the broken phase analysis. However, the determination of the coupling in the symmetric 
phase had a large statistical error. This problem has been eliminated by one of us [15] who 
carried out a complete analysis modeled on the work of [3]. The numbers of [5] survive 
this test reasonably. Of course, the analytical method suffers from systematic errors of 
a different type. It is difficult to know how much of the discrepancies at larger mass 
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values in the broken phase between the analytical and the Monte Carlo results are due to 
either approach. Still, the comparison gives a worst-case scenario for the systematics in 
both cases, since the methods are so different that systematics are unlikely to conspire to 
work in the same direction. Indeed, while the discrepancies on the F4 lattice and on the 
hypercubic lattice [6] are roughly of similar magnitude, the sign is opposite. We cannot 
be very precise, but the general conclusion we draw is that for ^j- as a function of mjj 
the overall systematic error in the Monte Carlo results cannot exceed a few percent for the 
highest Higgs masses measured. 

For £2 and S3 an analysis similar to [15] would be very demanding and has not been 
carried out. In [5] and for the lighter Higgs masses of actions S2 and £3 the Higgs mass was 
extracted using only the a field correlation function. Since now, for the larger masses, the 
width is suspected to be larger some additional checks are needed. The most direct test is 
to write down reasonable operators that would create predominantly two pion states and 
check for mixing effects. We did this for action S\ at /3q = 0.10, the coupling from which 
the bound quoted in [5] was extracted, and also for actions £2 and £3 at the couplings 
where we extract the bound. The large N results of section 2.3 show that more ambitious 
attempts to actually measure the width on the lattice by making the pions massive [11] 
may not be directly relevant to the massless case. 

The mass is now estimated by measuring the correlation matrix C(t) 

Cij(t) =< OiftOjiO) >-< Oi(t) >< Oj(0) > 

between several operators, 0$, described in the previous subsection, evaluated on two time 
slices separated by t lattice spacings. The eigenvalues are determined by diagonalizing 
C(to) -1 ' 2 C(t) C(to) -1 / 2 as in [16] with to = 1 except for the largest lattice with action 
Si where to = was used. After some trials it was decided to diagonalize a 3 x 3 matrix. 
For each lattice size we evaluate the lowest eigenvalue with this method. We also compute 
a "trial mass" from the a field correlation function alone. When the two numbers are 
different within errors we expect to have some level repulsion between the two lowest 
eigenvalues. We then identify the eigenvalue closer to the free two pion energy as the 
energy of the two pion state and the other as the resonance energy. To approximately 
correct for the repulsion we adjust the resonance energy by the amount that the two pion 
energy differs from its free value. This technique should eliminate leading 1/L corrections 
when they are significant numerically. The so obtained Higgs masses are extrapolated 
to L = 00 by assuming an behavior. The quality of these fits is acceptable but the 
motivation for the method of analysis is somewhat empirical and therefore we allow for an 
order 3% systematic error in the mass determination. 

The low lying spectrum for selected high mass values for each action are shown in 
Figures 3.4, 3.5 and 3.6 with the actual numbers given in Tables 6, 7 and 8. The new 
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Figure 3.4 The low lying spectrum for action S\ for various lattice sizes as measured in 
a numerical simulation at (3q = 0.10. The dotted lines correspond to the two 
lowest energies of two free pions. 

method of analysis confirms that the older results for action S\ are acceptable, which in 
turn increases our confidence when the method is applied to actions £2 and S3. We see 
no indications for strong mixing because the largest level repulsions observed are small 
compared to the eigenvalues themselves. The two-pion states have a spectrum well de- 
scribed by the lattice free particle dispersion and show only small amounts of sensitivity 
to the resonance. Taking into account statistical errors we settle on an error estimate on 
the mass determinations of 5%. This results in an error of about 6% on the determination 
of the physical Mfj. 



4. MAIN RESULT. 



A realistic and not overly conservative value for the Higgs mass triviality bound is 
710 GeV. At 710 GeV the Higgs particle is expected to have a width of about 210 GeV 
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Figure 3.5 Same as Figure 3.4 but for action S2 at (3q = 0.12, /?2 = —0.11. 



and is therefore quite strongly interacting. We estimate the overall accuracy of our bound 
as ±60 GeV. This includes the statistical error and systematic uncertainty of the measure- 
ments of Mjj (see section 3.4) as well as the systematic uncertainty assigned to fine tuning. 
The latter should also allow for effects of the hereto neglected 62 coupling in eq. (2.1.1). 
The meaning of the error is that it would be quite surprising if evidence were produced 
for a reasonable generic model of the scalar sector with observable cutoff effects bound by 
a few percent, but a Higgs mass larger than 770 GeV. It would be even more surprising if 
a future analysis ended up concluding that the bound is some number less than 650 GeV. 

The first estimates for the nearest-neighbor hypercubic actions gave a bound of 640 GeV 
[3,4] which is about 10% below the new number whereas the old F4 bound was 590 GeV 
[5]. Thus, the older results have proven to be quite robust and this is an indication that 
more search in the space of actions is unlikely to yield surprises. 
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Figure 3.6 Same as Figure 3.4 but for action £3 at (3q = 0.1175, = —0.11. 
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APPENDIX A. 



For the interested reader we sketch in this appendix the large N analysis for action S3 
that was not treated in [9]. We can write the action S'^ as in eq. (5.1) of [9] 



^ Jx,y ° e Jx 



(A.l) 



with 



9x,y = Q$x, y - - E + 12 E ^VJ 

x' n.n. to x x" n.n.n. to x 

g(p) = o Yl I 2 ~ cos (^m + vv) - cos(p M - p v )] 



(A.2) 



~Y2 E[ 2 - 2cos ( 2 ^)]+ E 
I » {^=±1} 

P 2 + 0( P 6 ) . 



1 - cos ( E 



Here we used S x ,y = l/2d Xj y for an F4 lattice and B* is its Brillouin zone. 

With the above notations the computations of the phase diagram, Higgs mass and 
cutoff effects are just as in sections 5 and 6 of [9]. We only need a few numerical constants. 
They are defined in [9] and for the present case take the values 

r = 0.09932603 

ci = 0.0281844 , c 2 = 2.2639 • 10~ 4 
7 = 0.01089861 
C=0 • 



(A.3) 



The resulting large N plots for the Higgs mass and cutoff effects, analogous to those in [9] 
for the actions considered there, are shown in Figure A.l. 



APPENDIX B. 



In this appendix we collect the numbers obtained from our simulations in several tables 
matching the graphs shown in the main text. 
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Figure A.l (a) mjj/ 'f n vs. m#: The solid line corresponds to = an d the dotted 
line to the optimal value = — @2,t.c- ( D ) Leading order cutoff effect in the 
width to mass ratio, (c) Leading order cutoff effects in the invariant tv — tv 
scattering amplitude at 90^. Here the dotted line represents center of mass 
energy W = 2Mjj, the dashed line W = 3M# and the solid line W = AM}j. 
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Table 1 



02 


A)c 


0.000 


0.0917( 2) 


-0.040 


0.0995( 5) 


-0.080 


0.1070(10) 


-0.110 


0.1118( 3) 


its of the critical line for action S2- 


02 


Poc 


0.000 


0.0653(3) 


-0.040 


0.0825(5) 


-0.080 


0.0990(4) 


-0.110 


0.1113(2) 



Table 2 Some points of the critical line for action S3. 



A) 


m H 


M H (GeV) 


0.0925 


0.198(37) 


438(84) 


0.0950 


0.414(25) 


524(32) 


0.0975 


0.550(27) 


546(27) 


0.1000 


0.702(30) 


593(25) 


0.1050 


0.830(41) 


576(30) 


0.1100 


0.946(50) 


568(30) 



Table 3 The Higgs mass rrifj in lattice units and the Higgs mass Mjj = x 246 GeV 
in physical units from the numerical simulations of action S\. The errors 
quoted are statistical errors only. 



A) 


m H 


M H (GeV) 


0.1120 


0.126(15) 


529(62) 


0.1130 


0.239(13) 


571(30) 


0.1150 


0.391( 8) 


627(12) 


0.1175 


0.534( 8) 


672(10) 


0.1200 


0.646( 9) 


696(10) 


0.1225 


0.721(10) 


694(10) 



Table 4 Same as in Table 3 but for the action ^2 at = —0.11. 



A) 


m H 


M H (GeV) 


0.1125 
0.1150 
0.1175 


0.302( 7) 
0.528(12) 
0.663(11) 


634(15) 
699(16) 
708(12) 
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Table 5 Same as in Table 3 but for the action S3 at fa = —0.11. 



L 


Three lowest energy levels 


Two lowest free n — n 


8 


0.802(24) 


1.202(54) 


1.818(77) 


1.362 


1.571 


10 


0.801(25) 


1.005(47) 


1.315(67) 


1.089 


1.257 


12 


0.754(13) 


0.894(18) 


1.063(19) 


0.907 


1.047 


14 


0.739(16) 


0.804(23) 


0.916(27) 


0.777 


0.898 


16 


0.721(33) 


0.679(23) 


0.839(23) 


0.680 


0.785 



Table 6 The low lying spectrum for action S\ for various lattice sizes as measured in 
a numerical simulation at fa = 0.10. The two lowest energy levels of a free 
two pion state of zero total three-momentum are also given. 



L 


Three lowest energy levels 


Two lowest free n — n 


8 


0.792(24) 


1.277(102) 


1.711(223) 


1.362 


1.571 


10 


0.775(14) 


1.115( 36) 


1.171( 39) 


1.089 


1.257 


12 


0.686(16) 


0.889( 32) 


1.034( 38) 


0.907 


1.047 


14 


0.673(15) 


0.794( 24) 


0.914( 25) 


0.777 


0.898 


16 


0.656(20) 


0.728( 21) 


0.863( 27) 


0.680 


0.785 



Table 7 Same as in Table 6 but for action S2 at fa = 0.12, fa = —0.11. 



L 


Three lowest energy levels 


Two lowest free n — n 


8 


0.817(15) 


1.250(146) 


1.660(196) 


1.369 


1.571 


10 


0.759( 9) 


1.108( 53) 


1.267( 64) 


1.091 


1.257 


12 


0.740( 8) 


0.964( 16) 


1.050( 15) 


0.908 


1.047 


14 


0.684( 8) 


0.817( 14) 


0.913( 15) 


0.778 


0.898 


16 


0.659( 8) 


0.719( 10) 


0.811( 11) 


0.680 


0.785 



Table 8 Same as in Table 6 but for action S 3 at fa = 0.1175, fa = -0.11. 
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